Quantum simulation of the Klein paradox with trapped ions 
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We report on quantum simulations of relativistic scattering dynamics using trapped ions. The 
simulated state of a scattering particle is encoded in both the electronic and vibrational state of an 
ion, representing the discrete and continuous components of relativistic wave functions. Multiple 
laser fields and an auxiliary ion simulate the dynamics generated by the Dirac equation in the 
presence of a scattering potential. Measurement and reconstruction of the particle wave packet 
enables a frame-by-frame visualization of the scattering processes. By precisely engineering a range 
of external potentials we are able to simulate text book relativistic scattering experiments and study 
Klein tunneling in an analogue quantum simulator. We describe extensions to solve problems that 
are beyond current classical computing capabilities. 
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Simulating quantum mechanics using conventional 
computers rapidly becomes intractable as the physical 
systems to which it is applied grow larger. A proposed 
solution is to use highly controlled laboratory quantum 
systems themselves to perform such simulations [TH3]. 
Systems currently under investigation for this purpose 
include photons 01 [5], trapped atoms [6l [7] and su- 
perconductors [8]. A particularly promising approach 
uses trapped ions [5Hll) with which several quantum 
simulations have recently been performed. One line of 
work is the simulation of quantum models of interacting 
spins [ini [12 [T3]. Here, the internal states of ions en- 
code the spin states and spin-spin interactions are simu- 
lated by laser induced state-dependent forces. Recently, 
a quantum simulation of the dynamics of a free rela- 
tivistic particle has been performed in our group |14j . 
In this case, in contrast to simulations of spin systems, 
both discrete and continuous variables have to be simu- 
lated. In this paper we perform quantum simulations of 
the scattering dynamics of relativistic quantum particles. 
We show how a wide range of external potentials, from 
which the simulated particle can scatter, is engineered 
using two ions coupled via laser fields. We measure wave 
packets of the scattering particles, visualizing Klein tun- 
neling 'frame-by- frame'. Finally we describe extensions 
to efficiently simulate processes that are beyond current 
classical computing capabilities. 

In its original form [15] the Klein paradox considers a 
relativistic electron described by the Dirac equation, with 
total energy E and rest mass energy mc^, hitting a step- 
shaped potential barrier of height V . For barrier heights 
smaller than the kinetic energy of the electron, V < 
E — mc^, the particle is predicted to partially transmit. 
For a slightly larger barrier, E + mc^ > V > E~mc^, the 
particle should completely reflect. These situations agree 




FIG. 1: (a) Klein tunneling in position i) and momentum 
space ii), for a relativistic particle scattering from a linear 
potential (j>{x) = 9^- A wavepacket initially in the positive 
energy branch and with positive momentum (I) moves up a 
slope where it can reflect (II) while remaining in the positive 
energy branch or tunnel (III) while switching energy branch. 
(b) The Dirac spinor is encoded in an internal state of ion 1 
and a collective motional state of both ions. The Hamiltonian 
for a free Dirac particle {Hd) is implemented by bichromatic 
laser 1 which couples ion 1 to the collective motional mode. A 
second bichromatic laser couples ion 2 to the same motional 
mode and in this way creates the potential (j>{x) (color online) . 

with the predictions of the (non-relativistic) Schrodinger 
equation. However, for V > E ^ mc^ the particle can 
propagate undamped in the potential barrier, by turn- 
ing into its anti-particle. This effect is known as the 
Klein paradox. Klein's results have been extended to 
other types of potentials [TMTB] and the physics of Klein 
tunneling emerges for electrons in graphene [19H21| . In 
quantum field theory, the paradox is resolved by the no- 
tion of pair creation by the external potential [551 [23] ■ 



In one dimension, the Dirac equation for a particle in 
an electrostatic potential 4>{x) is given bv|24) 



ih 



dt 



[cpax + rnc a^ 



(x)h) *. 
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Here c is the speed of light, p the momentum opera- 
tor, m the particle mass and e its charge. The matri- 
ces (Tj, j — x,y,z are the Pauli matrices and I2 is the 
identity matrix. In one dimension there is no spin (no 
magnetic fields) and therefore the wave functions ^ are 
2-component spinors, reflecting that there are positive 
and negative energy solutions E± — zty'c^p'^ + ra^d^. A 
spinor allows for arbitrary superpositions of these com- 
ponents. Free particles {4>{x) = 0) in one of the two en- 
ergy branches remain there indefinitely, but for 0(a;) 7^ 
the spinor can switch energy branch and Klein tunneling 
can occur. For linear potentials e(j){x) — gx the situa- 
tion is conceptually equivalent to Landau Zener tunnel- 
ing il6j ,24j .25, , as shown in FigfTJa). Depending on the 
size of the splitting of the two energy branches (2toc^) 
and the acceleration (17/771), the particle can either adi- 
abatically follow the positive energy branch and reflect, 
or make a non-adiabatic transition to the negative en- 
ergy branch and tunnel. For ultra-relativistic particles, 
the probability for tunneling is given by Ptunnei = e"^'^'", 
with r = m?c^ /2hg. 

In our experiment, two ^'^Ca^ ions are trapped in 
a linear Paul trap with trapping frequencies LOax — 
2tt 1.36 MHz axially, and uJrad — in 3 MHz radially. 
A spinor is encoded by mapping the continuous posi- 
tion and momentum components to the quadratures of 
a vibrational mode, and the discrete components to two 
long-lived internal electronic states of ion 1 (Fig. flVb)). 
Evolution under the free Dirac Hamiltonian is realised 
using a bichromatic light field [TH [5^ with an over- 
all detuning coupling the electronic state of ion 1 to 
the vibrational mode via both red and blue sidebands. 
Within the Lamb-Dicke approximation, this interaction 
yields the Hamiltonian Hd = 2T]Afliax p + hfliaz , 
where rj — 0.044 is the Lamb-Dicke parameter, fii is 
the bichromatic Rabi frequency and 2fli is the overall 
detuning. The momentum operator p = 7/7(0^— a)/2A, 
where A = ^yh/AmtUax ~ 7 nm is the size of the ion's 
ground-state wave function, 777 the ion's mass (not to be 
confused with the mass of the simulated particle) and a^ 
(a) are the creation (annihilation) operator on the centre- 

(n) 

of-mass vibrational mode. The notation cr) represents 
Pauli operators on the subspace of ion n. The mapping 
to the free Dirac Hamiltonian is complete by making the 
identifications c := 277J71A and mc^ := hili. 

External potentials are simulated using a second ion 
(ion 2) and another bichromatic light field coupling 
ion 2 to the same vibrational mode but via a differ- 
ent electronic transition (Fig. fTlb)). With appropri- 
ately set phases the interaction Hamiltonian is given by 



e0i = gai: x with g = h7]^2l ^ |21]- Here, D.2 is the 
Rabi frequency of the second bichromatic light field and 
x ^ {a) + a)A is the position operator. When ion 2 is 

prepared in an eigenstate of cri this operator can be re- 
placed by its -|-1 eigenvalue, reducing the interaction to 
a linear potential and F — ^\/{ArfQ.iQ,2)- 

To reconstruct the spatial probability distribution of 
a spinor after the simulation, another bichromatic light 
pulse maps information about the position/momentum 
of the spinor onto the internal state of ion 2 [T3J [?7] . 
A measurement of this state, as a function of the pulse 
length, can be used to reconstruct the probability dis- 
tribution (see reference [28] and its online material). In 
short, (1) after the simulation ion 2 is prepared in one 
of the states ^(J), or (1)2' (^) ^ bichromatic light field 
is used for a time tprobe to implement a displacement 
operation of the form Uprobe — exp(— j/cxcTy /2), with 
k ex tprobe- This Operation causes the motional state to 

split in phase space along the momentum- axis. The re- 

(2) 
suiting interference fringes in the observable ai , which 

can be measured by fiuorescence detection, as a func- 
tion of tprobe then represent the fourier components of 
the original motional state sin(fca;) or cos(A:x), depending 
on which state ion 2 was prepared in. A fourier trans- 
form (3) of this data gives access to the spatial probabil- 
ity distribution. Before this measurement, the internal 
states of the two ions can be traced out by a series of 
laser pulses pumping the internal states into (j^) and 
{i)„- For spinors with high momentum, it is also pos- 
sible to obtain the wavepackets associated with positive 
and negative energy separately. In these cases the in- 
ternal spinor states are not entangled with the motional 
state, and are given by the states (_|^^) . A 7r/2 pulse can 
be used to map either of these energy states to (J) ^ after 
which a short (200 /is) projective fiuorescence measure- 
ment is done. Since the state (J)^ scatters no photons 
during fiuorescence detection, the motional state of the 
ions remains intact and can be analyzed afterwards. The 
cases where ion 1 is found in (°)^ are discarded. 

At the start of each experiment, laser cooling, opti- 
cal pumping and coherent laser pulses, in a magnetic 
field of 6 G, prepare the ions in the axial center-of- 
mass mode ground state and internal states 

and ^fi)„ (see Fig. 1.). We create an initial spinor 



\/2Vl)l 

75(1)2 (see Fig. ] 
state corresponding to a particle with (x) = and 
Po — iP) — 3.5/7/A by a displacement along the mo- 
mentum quadrature. The generated spinor, 'ijj{x; t = 

0) ex e^P°^/'^e~i^ (J)^, is comprised largely of positive 
energy: \{ip{x;t = 0)\P+\'iP{x;t = 0))\^ > 0.98, where 
P+ is the projector onto the positive energy state. The 
rest mass energy is set to hfli — h2TT 1.3 kHz and fli — 
27r 17.5 kHz corresponding to an equivalent speed of light 
of c ~ 0.01A//is. Fig. [2]shows the scattering for different 
slope gradients (g), achieved by setting f22/27r = 0, 22, 50 




X{^) " l(p)l(ft/A) 

FIG. 2: Quantum simulation of relativistic scattering for linear potentials. Measured particle wave packets (filled curves) 
are compared with ideal predictions (solid black lines) and predictions taking corrections to the Lamb-Dicke approximation 
into account (dashed black lines). In the first and last frames of each sequence the positive (green) and negative (red) energy 
components are reconstructed separately. The blue color scale of these panels represents the measured expectation value of 
momentum. The axis on the right shows the potential energy in units of initial kinetic energy. Without a potential, the particle 
moves to the right with constant velocity (a). For a shallow potential gradient, the particle is almost completely reflected (b) 
and for steeper gradients (c,d) part of the wave packet propagates into the repulsive potential via Klein tunneling (color online) . 



and 76 kHz, corresponding to Ptunnei — 0, 0.03, 0.21 and 
0.36 respectively. The measured tunnel probabilities are 
given by 0.017(7) {0}, 0.10(1) {0.07}, 0.32(2) {0.22} and 
0.45(3) {0.39}, here we put results obtained by numerical 
calculations in curly brackets. Animated versions of each 
of the cases can be found in the online material linked 
to the paper (EPAPS Document No. [number will he 
inserted by publisher]). 

We are also able to create approximately quadratic 
electric potentials which, due to Klein tunneling, are non- 
confining in the relativistic limit [18J. The quadratic po- 
tential is implemented experimentally by detuning the 
bichromatic beam on ion 2 by 21^2 1 so that the cou- 
pling between this ion and the vibrational mode becomes 



The scattering dynamics in this case correspond to that 
of a quantum relativistic mass on a spring that is given a 
small initial kick and show Klein tunneling at each turn- 
ing point. 

The state used for the simulations in Figs. l3[a) & 
(b) was created by applying a bichromatic displacement 



operation to ion 1 of the form Up 



-iH„ 



^t/h 



with Hprep2 — hr]fl 



.(1), 



.(2); 



J2) 



hr}Q,20'x x/A + hfl2a-z ■ In the limit of a large detuning, 
^2 ^ ??^2 the effective interaction Hamiltonian becomes 
with q=hr]'^nl/ (2A'^n2) I24j. Preparing 



e(P2 



(2) -2 



ion 2 in the -1-1 eigenstate of cr^ this reduces to a quadratic 
electric potential. 

Figure |3] shows results from our investigation of 
quadratic potentials. The first sequence (a) shows a 
particle, with initially positive energy, {x) = and 
(p) = 0, evolving as a free particle {^2 = 0). As ex- 
pected, the wave packet simply disperses. Sequence (b) 
shows the same initial state evolving under a quadratic 
potential generated by setting Q.2 = 27r 33 kHz and 
^2 = 27r 50 kHz, such that 77^^2/^22 = 0.067 < 1 and 
q = 27r 73 Hz hj t^ . The dynamics are clearly different 
in this case: the wavepacket is still unconfined, due to 
Klein tunneling, but it spreads more slowly. Sequence (c) 
shows results for the same potential, but for a positive en- 
ergy particle with an initial momentum (p) = 0.23 /i/A. 



^prep2 — ^ 

prep^^o'x ' / A. and ion 1 prepared in 
the internal state (0)2- ^ 16 /iS pulse with VLprep2 — 
2tt 83 kHz produces a spinor with (p) — and \{'ip{x;t — 
0)\P+\^j{x;t = 0))|2 > 0.98. The state used in Fig.|3](c) 
was created in a similar way, but an additional pulse was 
used to give the state a momentum of (p) = 0.23 h/A. 

There are a number of errors that reduce the quality of 
our quantum simulations. The internal state coherence 
time for the ions is ~ 3 ms, limited by magnetic field fluc- 
tuations, while the whole simulation takes up to 1.5 ms. 
The motional coherence is limited by slow drifts, which 
can change the trap frequency by about « 25 Hz. Both 
effects cause broadening of the wave packets, whereas sys- 
tematic errors in state preparation due to slowly varying 
experimental settings can cause additional structure in 
the reconstructed data. For the steepest slope, states 
of more than 150 phonons on average were created, for 
which the Lamb-Dicke approximation starts to fail. This 
changes the simulated Hamiltonian somewhat, while also 
affecting the reconstruction [25]. All these errors can be 
improved by technological development. 

Scaling up our quantum simulation is conceptually 
straightforward: additional Dirac particles can be sim- 
ulated by adding more ions and motional modes. Each 
extra spatial dimension requires an additional motional 
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FIG. 3: Quantum simulation of scattering for quadratic 
potentials. Measured particle wave packets (filled curves) 
are compared with theoretical predictions (black lines) for 
quadratic potentials. The particle's rest mass energy is set to 
h£l\ — h2TT 0.65 kHz. In (a) & (b) the initial state has almost 
purely positive energy, and momentum (p) — Oh/ A. Without 
a potential (a) the wave packet diffuses. With the potential 
switched on (b) the spreading is significantly restricted, but 
can continue through Klein tunneling. In (c) the initial state 
has a small momentum (p) = 0.23ft/ A and the wavepacket os- 
cillates. The red dots represent (x). The lower three figures 
show larger versions of some of the wave packets. The axis 
on the right is in the same units as in Fig. 2. (color online) 

mode. Performing classical simulations at the same reso- 
lution achieved by the quantum simulations requires de- 
scribing a 2 X 200 dimensional Hilbert space (one qubit 
and 200 harmonic oscillator states). Adding just one 
Dirac particle increases the Hilbert space to ~ IGxlO** 
dimensions, which would require 10'^ GB of memory just 
to store the Hamiltonian operator with 64-bit precision. 
Solving the scattering dynamics of two interacting Dirac 
particles is in general an open problem into which an ion 
trap quantum simulator could provide new insight. It is 
also possible to extend this work to simulate interactions 
with magnetic, scalar and much studied confining po- 
tentials such as the Dirac oscillator [35]. Other possible 
extensions would be to simulate quantized Dirac fields or 
Majorana physics pO] . 
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